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We discuss a novel form of the variational approach in Quantum Field Theory in which 
the trial quantum configuration is represented directly in terms of relevant expecta- 
tion values rather than, e.g., increasingly complicated structure from Fock space. The 
quantum algebra imposes constraints on such expectation values so that the variational 
problem is formulated here as an optimization under constraints. As an example of ap- 
plication of such approach we consider the study of ground state and critical properties 
in a variant of nonlinear sigma model. 

The variational approach is one of the corner-stones of nonperturbative methods 
in Quantum Mechanics and Quantum Field Theory. In this approach, the expec- 
tation value of the Hamiltonian is analyzed on a set of quantum configurations of 
specific form and its minimum is sought. Variational Method takes roots in Ritz 
Theorem [1] which states that for a hermitian Hamiltonian operator with the spec- 
trum bounded from below 

($\H\$) 

<"> - w - E °- (1) 

for any quantum state Here Eq is the lowest eigenvalue of H . Ritz theorem can 
be transparently motivated using the fact that hermitian operator should have a 
complete set of eigenstates, i.e. an arbitrary quantum state can be represented as 
a linear superposition of the Hamiltonian eigenstates. Then, after doing a simple 
algebra one can get 

=X>nMVfc|tfhM =Y J \'t>n?E n {^ n ) >E (<i>\<P}. (2) 

n,n' n 

Here |$) = ^2<f>n\ipn) an d \ip n ) are eigenstates of Hamiltonian H with eigenvalues 

n 

E n > Eq, From this simple argument, it follows that one will get an upper limit 
for the ground state energy if one computes the Hamiltonian expectation value on 
arbitrary quantum state $. This suggests a way to generate a useful estimate for 
the ground state energy by considering a class of trial quantum states {<& a , ol G A} 
and minimizing E(a) = ^"^^"^ with respect to parameter a. 

The variational approach found a wide range of applications due to its simplic- 
ity, possibility of analytical results and important nonperturbative content. One of 
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the simplest examples is the problem of Hydrogen-like atom which is described in 
dimensionless unitsby the Hamiltonian 

1 9 1 

ff=--V 2 --. (3) 

Choosing an exponential trial wavefunction $ Q (r) = e~ Qr , one trivially obtains 

{H) = \a{a-2) (4) 

from which the solution a = 1, E = —0.5 immediately follows. In this case, the 
variational estimate is in fact the exact answer. 

Many examples of the original applications of variational method in Quantum 
Mechanics for computation of the energy levels of simple molecules can be found 
in the literature. In the later time, variational method had been employed increas- 
ingly for the precision many-body calculations of the ground state energies and 
first excited levels of atoms and simple molecules. In theoretical studies, varia- 
tional method is often used to derive important nonperturbative constraints, e.g. 
asymptotic properties, on quantum-mechanical models. First work using variational 
method in Quantum Field Theory dates back to 1960th since it had been applied 
to a variety of atomic, nuclear and quark problems. In QFT, variational method is 
commonly used to investigate stability of the vacuum as well as to study dynam- 
ics of one- or two-particle excitations [2] . Also, the conclusion about possibility of 
color-superconductivity in QCD vacuum can be justified by considering the trial 
states with quark-quark vacuum condensations. The use of variational principle to 
formulate nonperturbative equations of state for few-body systems had also enjoyed 
growing popularity [3,4]. 

In QFT, the variational principle is derived from the equation for the covariant 



= M 2 |$), (5) 

where P^ is the momentum-energy operator. In the Center of Momentum Frame, 
P|<6) = 0, this becomes 



(P°) 2 |$) = M 2 \<S> 
p°m = ±M\$). 



(6) 



Eq.@ describes two energy spectra identical up to the change of sign which are 
conventionally associated with particle and antiparticle sectors of a field-theoretical 
model. By excluding spectrum of negative energies, one arrives at 

or, equivalently, 



(5<f>\H - = 0. (8) 
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Variational principle JSJ applied to the states from many-body Fock space 

l*>=E /// dki...dk n C (i) (k 1 ,...,k n )a{ i ...a{jO) (9) 



typically leads to an infinite system of coupled equations for the amplitudes 
C7j)(fei, . . . , kn). This is usually simplified using Tamm-Dancoff truncation [5] or 
a variational ansatz [3]. Ansatzes involving a handful of few-body excitations, e.g. 

I*) = (1 + / dkf k a\ +JJ dhdhfk^alaljO), 

are often used in addition to ansatzes defined by coherent states 

|$) = cxp( / dkf k al)\0) (10) 



|$) = exp( J f dhdk 2 f klk2 alal 2 )\0). (11) 

While variational approach in QFT retains such important features as simplicity 
in application and formulation as well as significant nonperturbative contents, it 
suffers from numerous drawbacks which hampered its applications in QFT. Among 
these is the lack of explicit covariance, the necessity of nonperturbative regulariza- 
tion/renormalization if the model is not a-priory finite and the growing complexity 
of trial quantum states. In this paper, we are going to consider a new variational 
approach that may be helpful in dealing with the latter problem. 

We note that, while in Quantum Mechanics the quantum state could be de- 
scribed by a single wavefunction, in QFT one needs an infinite number of many- 
body amplitudes CuUki, . . . , k n ), as introduced in Eq.(PJ), to properly describe the 
quantum state. This makes general variational problem intractable and even simpler 
variational ansatzes, like in Eas. (|10|l and (jl 1|> . require significant algebraic work to 
obtain higher particle number contributions. In QFT, one has to deal with an infi- 
nite number of many-body amplitudes even though one may only need to compute 
a limited number of the relevant expectation values, e.g. (tt 2 ), ((V</>) 2 ), (</> 2 ), (<p 4 ) 
in </> 4 scalar field theory. This overcomplicated nature of variational approach in 
QFT led us to the following idea. Instead of going through the complex Fock space 
representation of the trial quantum state in computing the expectation values of the 
relevant operators, one may parametrize the trial configuration in terms of the ex- 
pectation values themselves, e.g. = (tt 2 ), A^2 = (<fi 2 ) and A^4 = (0 4 ). Because 
of the quantum nature of operators, certain constraints would be imposed on the 
expectation values A's by quantum algebra. Then, the variational problem is cast 
into a constrained minimization problem in terms of only the relevant expectation 
values, e.g. 



+ A V + m 2 A 02 ) + Aa 04 - 
constraints on A„, Av^, A^2, A^4 
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This idea is specifically motivated by the following theorem. 
First Symmetric Decomposition Theorem: Symmetrized Fock space T = 
{| 77)} is isomorphic via 

(v\a{a k '\v} =9{k,k') 

to the space of linear integral operators with kernel g(k, k') that are positive definite 
hermitian and have finite traces. By this we claim that, if g(k, k') is known to define 
a hermitian positive definite linear integral operator with a finite trace, then there 
exists a normalized quantum state in symmetrized Fock space T which gives g(k, k') 
as (a^,afc/) and vice versa. 

To prove the First Symmetric Decomposition Theorem, we note that for hermi- 
tian operator introduced by kernel g(k,k') there should exist a full orthonormal set 
of eigenfunctions and eigenvalues such that 

g(k,k') = J2^m9m(k)g* m (k'), (12) 

m 

where g^ik') is a complex conjugation of eigenfunction g m {k) and for positive def- 
inite operators the eigenvalues A m > 0. It is sufficient to consider a set of states 

\n) = ^=Jjdk 1 ...dk n f*(k 1 ,...,k n )a{ i ...a{ n \0), (13) 

where 

fn(ki, ...,k n ) = ^2 \/^9m( fc i) • • • 9m{k n )- (14) 

Given orthogonality and normalization of the eigenfunctions, it is easy to check that 
{n\a\a k , \n) = ^ ^m9m{k)g* m {k') = g(k, k'). 

m 

Then, it is always possible to choose a set of amplitudes fi n such that for \rj) = 

n 

{■q\a\a k ,\ri) = g(k,k') and (r]\r]) = 1. 

Note that a finite trace is required provided that 

( V \N\ V )=Ti[g(k,k')]<oo. 

It is trivial to show that for any \r]) from the Fock space the expectation value 
{a\dk') defines a positive definite hermitian linear integral operator with a finite 
trace. This completes the proof. 

Furthermore, it is possible to extend this statement to Second Symmetric 
Decomposition Theorem stating that for expectation values on a normalized 
quantum state from symmetrized Fock space |?7), 

A v (k, k') = {r)\a\a k > + tf_ k ,b- k \r)), 
B v (k,k') = {i 1 \a\bl kl +a k ,b_ k \r l ) 1 
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treated as kernels of hermitian linear integral operators, precisely 



1 + A v y Jl + B%, (16) 



i.e. 1 + A n — 1 + B 2 is a positive definite operator a . 

The significance of these statements is in the claim that any expectation value 
of the form given above can be represented in terms of a state from the Fock space 
and vice versa. Thus, in dealing with these expectation values one need not an 
explicit complex structure from the Fock space but may work entirely in terms of 
the relevant expectation values constrained by a condition like in Ea. H16|) . 

We will now illustrate application of this principle to the example of a study of 
ground state and critical phenomena in a variant of nonlinear sigma model. Advan- 
tage of this model for our application will be that it is completely formulated in 
terms of the expectation values of the quadratic operators in and, thus, our re- 
sults about isomorphism between T and (A rp B n ) can be straightforwardly applied. 
O(N) nonlinear sigma model is introduced as a free field theory constrained to live 
on a sphere of radius R. It is defined by a free field Hamiltonian 

H = [ dx±(f + + if?) (17) 



2 

and a constraint for ./V-component vector <j>(x) 

\$(x)\ 2 = R 2 . (18) 

In our model we will enforce this constraint softly on average, so that in the canonical 
quantization Eas. i|17fl and i|18|) will be given by 

(:ff:)~J& fc E4(M) 

1 frK (19) 

(: |0 2 (x)| :) ~ 5dKdk- J §—Y J [A\{KK + k) + Bl 1 {k,K + k)] 1 ' 

i 

We define the physical subspace as the subset of the Fock space T satisfying condi- 
tion {t\\ : \c/> 2 (x)\ : — R 2 . Consequently, we define the ground state as the state 
from the physical subspace with the lowest energy. 

We reformulate an original variational problem for the ground state in terms of 
the expectation values themselves, subject to the quantum algebra constraint Ijl6|l . 
In these terms for the hermitian linear integral operators and B v , we have 

Tr[e • A % S\ — > min 

Tr[e- 1 -(Aj,+4)]=i2 2 (20) 
off-diagonal / dk-^= EI^O' K + k ) + B i,(k, K + k)] = 0, K ^ 0. 



Here e is the diagonal matrix with = yjk 2 + ^ 2 5kk' , " ■" stands for matrix 
multiplication and Tr stands for integration over momenta and summation over i. 

a notc that = (r)\B\rj) {ri\B\rj) ^ {v\B 2 \rj) for B = aW + ab 
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From Ea. (|16fl . we know that for given B^ linear operator A 1 can be represented by 



(21) 



A; = M l + yjl + (B^) 2 - 1, where M l t 0. Then, 
Tr[e • At] — » min <^ 



Tr[e • M l ] + Tr[e ■ (^1 + (fi*) 2 - 1)] -> min . 



Given M. > 0, diagonal elements M.\ k > and Tr[e • M 1 ] > so that to minimize 
the expectation value of the Hamiltonian one needs obviously M. % — > 0. Thus, the 
constraint Ea. (|16|l is resolved and we get 



Tr[e- (^l +^- l)]^ min (a) 

Tr^ 1 • (^1 + (.5 ^-1 + 4 )] = R 2 (b) ( 22 ) 



One can observe that for the solution of Eas. l|22h ,) and (|22b ) B^ should be diagonal. 
For off-diagonal k ^ k' , one may consider variation 5B kk , that preserves Eq.JSJb) 
and show that 




<5Tr[e • Ul + (Si) 2 - 1)] = (1 - -f )ey , " 6B\ k „ (23) 



which implies that B^ with off-diagonal nonzero elements cannot be the solution of 
the minimization problem Ij22(l . In that case Eq. I|22h ) is redundant as it is automat- 
ically satisfied by the solution of Eqs.(j22t) and 122b ). 

Having said that, the problem defined by Ens.l|22h,) and l|22b ) can be solved 
using the method of Lagrange multipliers. We introduce a new variable 



Ql,(k) = yfl + (B^k, k)) 2 - 1 + B;(fc, k), (24) 
and rewrite Ea. H22|) in the form 

[Q?,(fc)] 2 



. J dfc l+Q-(fc) efc mm 
EfdkQ^k^^R 2 , 

i 

which is rewritten with Lagrange multiplier A as 

^( k )EJdk(g^e k -XQ^k)e k 1 

and Eq. i|25|) is then given by 



(25) 



(26) 
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Ea. (|27|l is the gap-equation for X(R 2 ). In 2+1 dimensions, this integral can be 
computed exactly to yield 



Thus, we arrive at the following solution of our original problem. The ground 
state in our model, defined as the state from the Fock space T with lowest energy 
(: H :) and satisfying (: \4> 2 {x)\ :) = R 2 , is described for R 2 < R 2 C = by the 

expectation value (: 4>i(k)* <j>i{k) :) = , £ , fc — 1, where A is the solution of Ea.p8)l. 

V e 'k~ x 

For R 2 > R 2 the solution of the form (|26[1 can no longer be found because of 
the "finite capacity" of k ^ modes in the distribution (|26[l . Instead, the ground 
state is described as superposition of distribution l|26|l with A = fi 2 and a singular 
Bose condensation in k = mode. Development of Bose condensation for R 2 > R 2 C 
describes a phase transition of the second kind in this model. 

This discussion is the exact variational solution of the original problem thanks 
to the exact constraints on the expectation values and that we were able 
to find. As one can see, we rendered the explicit Fock space structure completely 
unnecessary. The final answer was given in terms of the expectation values of a 
given operator on the ground state without explicit reference to the Fock space. 

While the knowledge of an exact image of the Fock space T in terms of given ex- 
pectation values, say A v and B Vl is very interesting and, as we have shown, may be 
beneficial in certain problems, such detailed information for the expectation values 
of more complex operators, e.g. 4 , may be difficult to obtain. In general, if approx- 
imate constraints on the expectation values of the operators can be established, 
the approach outlined above would provide a variational estimate for the ground 
state of the model. Among its advantages would be elimination of the necessity to 
include the consideration of complicated states from the Fock space explicitly as 
well as a possibility to incorporate renormalizations via the expectation values of 
the quantum operators themselves as they enter the final constrained optimization 
problem. In this sense further investigations of this approach present clear interest. 
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